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Abstract Recently, ISun et all d2008l) published new Galactic 3D-models of magnetic 
fields in the disk and halo of the Milky Way and the distribution of cosmic -ray elec- 
tron density by taking into account the thermal electron density model NE2001 by 
ICordes & Laziol d2002l [2003b . The models successfully reproduce observed continuum 
and polarization all-sky maps and the distribution of rotation measures of extragalactic 
sources across the sky. However, the model parameters obtained for the Galactic halo, 
although reproducing the observations, seem physically unreasonable: the magnetic field 
needs to be significantly stronger in the Galactic halo than in the plane and the cosmic- 
ray distribution must be truncated at about 1 kpc to avoid excessive synchrotron emission 
from the halo. The reason for these unrealistic parameters was the low scale-height of 
the warm thermal gas of about 1 kpc adapted in the NE2001 model. However, this scale- 
height seemed well settled by numerou s investigations. Rece ntly, the scale-height of the 
warm gas in the Galaxy was revised bv lGaensler et al.l d2008l) to about 1.8 kpc, by show- 
ing that the 1 kpc scale-height results from a systematic bias in the analysis of pulsar data. 
This implies a higher thermal electron density in the Galactic halo, which in turn reduces 
the halo magnetic field strength to account for the observed rotation measures of extra- 
galactic so urces. We s lightly modified the NE2001 model for the new scale-height and 
revised the ISun et all I (120081) model parameters accordingly: the strength of the regular 
halo magnetic field is now 2 pG or lower, and the physically unrealistic cutoff in z for the 
cosmic -ray electron density is removed. The simulations based on the revised 3D-models 
reproduce all-sky observations as before. 
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1 INTRODUCTION 

It is important to obtain a realistic model of the Galactic magnetic field, the cosmic-ray electron and the 
thermal electron density distribution for the understanding of the physical processes in the magnetized 
interstellar medium. Relativistic electrons lose their energy in the Galactic magnetic field by polarized 
synchrotron emission. Synchrotron emission fluctuations are the primary contamination for the analysis 
of cosmic microwave background observations. The Galactic magnetic field causes deflections of ultra- 
high energy cosmic-rays, which must be corrected for to identify their origin. The magnetized thermal 
interstellar medium generates rotation measures (RMs) in the direction of extragalactic sources, which 
must be properly subtracted. However, it is very challenging to obtain realistic Galactic 3D-models be- 
cause of our position inside of the Galactic disk and the limited number of all-sky surveys needed for 
this task. 
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The results from earlier modelling o f the Galactic ma gnetic field, relativistic electron and thermal 
electron densities have been reviewed bv lSun etal 1 (12008b . Most of the models were derived from se- 
lected data sets only and do not agree with other observations. ISun et all d2008l) made the first attempt 
to establish 3D-emission models (SRWE08 models hereafter) aiming to properly represent all relevant 
radio observati ons available such as RMs of extraga l actic s ources, and total intensity and polarization 
all-sky surveys, [jansson et alj d2009h and lJaffe et alj d2010h foll owed that attempt an d additionally in- 
voked a quantitative comparison with the observati ons. However. IJansson et al. 1 J2009I) did not take into 
account total intensity maps and lJaffe et al.l d2010h so far presented a 2D disk model. ISun et al. I (120081) 
relied on a qualitative comparison between simulations and observations, which seems at the present 
stage sufficient to constrain 3D models describing the global Galactic properties. A quantitative com- 
parison needs to separate local large-scale features and more distant anomalies from the all-sky maps, 
which is a very ambitious task by its own and beyond our scope. 

The SRW E08 models are based on the thermal electron density model NE2001 bv lCordes & Laziol 
(120021 120031) . which uses a scale-height of about 1 kpc. To account for the RMs of extragalactic sources 
at high latitudes a strong regular halo magnetic field of up to about 10 fiG was required. Subsequently, 
the cosmic-ray electron density distribution needs to be cut in 43 at 1 kpc below and above the Galactic 
plane to avoid excessive polarized emission at high latitudes. However, both the strong halo field and 
the z-truncation of the cosmic-ray electron density are physically unrealistic. As already argued by 
ISun et all d2008l) this problem could be solved by increasing the thermal electr on density scale-heigh t 
by a factor of about 2. Recently, a scale-height of about 1.8 kpc was derived bv lGaensler et alj d2008l) . 
which is very close to our prediction. This motivates us to update the SRWE08 models. 

The paper is organized as follows: we briefly describe the method and available observations in 
Sect. 2, revise the SRWE08 models in Sect. 3, and summarize our results in Sect. 4. 



2 THE METHOD AND AVAILABLE OBSERVATIONS 

We follow the same method to obtain 3D-emission models as already detailed by ISun et "all d2008l) . 
The thermal electron density model N E2001 was the basis fo r the SRWE08 models and needs revision 
to adapt for the larger scale-height by iGaensler et al.l (120081) . The thermal electron density n e and the 
regular magnetic field component B« along the line-of-sight I were obtained from RMs of extragalactic 
sources, RM ~ n e Bnl. The random field and the cosmic-ray elect ron density were cons trained by 
total intensity and polarization all-sky maps. The Hammurabi code bv lWaelkens et al.l (120091) was used 
to simulate all-sky emission, which was then compared to the observations. Despite recent efforts the 
RMs of extragalactic sources are still sparsely distributed and are likely influenced by local structures 
on large scales in an unknown way, so that a quantitative x 2 -fit will not determine the halo parameters 
more reliably compared to our qualitative approach. 

We used RMs from the Canadian G alactic Plane Survey (CGPS, iBrown et all 120031) and the 
Southern Galactic Plane Survey (SGPS, IBrown et al.l 120071) to constrain the ma gnetic field in the 
Galactic disk. For the halo, previous RM measureme nts collected by Han et al.l (1 1997b and the Effels berg 



L-band RM survey (Han et al. in prep.) were used bv lSun et al. I (120081) . Recentlv. lTavlor et ail (12009b ob- 



tained RMs towards 37,543 extragalactic sources by re-processing the polariz ation data from th e NVSS 
survey dCondon et al.L fl998h . This is the largest RM data set available so far. iMao et al.l d2010b pointed 
out that individual NVSS RMs might have large errors, as they were derived from two frequencies only. 
The NVSS RMs are partly not reliable to wards the plane beca use the large RMs there are usually be- 
yond the ambiguity limits as explained by iTavlor et al.1 (12009b . However, this data set should describe 
Galactic large-scale structures well beyond the plane. Therefore we included these RM data to constrain 
the magnetic field in the halo. The RMs for the region of 100° < I < 120° and -5° < b < 20° from 

the CGPS extension (Brown et al. in prep.) were als o used. 

The 408 MHz total in tensity all-sky survey dHaslam et al.l 1 19821) and the WMAP five-years 



22.8 GHz polarization map dHinshaw et all 12009b both trace Galactic synchrotron emission and were 



1 The cylindrical coordinate (R, if>, z) is defined as: R is the Galactocentric radius, <f> is the azimuth angle starting from 1 = 0° 
and increasing in anticlockwise direction, and z is the distance to the Galactic plane. 
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used to constrain the total m agnetic field s and c osmic-ray electron density. The WMAP five-years MEM 
free-free emission template dGold et all 120091) served as a crosscheck for the thermal electron density 
distribution model. 

3 MODELLING REVISITED 
3.1 Thermal electron density 

NE2001 is a 3D-model describing the diffuse Galactic thermal gas dCordes & Lazidl2002l 120031 The 
thermal electron density distribution consists of several components: a thin disk, a thick disk, spiral 
arms and in addition a large number of thermal source complexes, which is, however, not complete. 
The ionized gas at high latitudes primarily attributes to the thick disk component with a scale-height of 
about 1 kpc and a mid-plane density of about 0.034 cm~ 3 . This scale-height was determined by fitting 
DM sin \b\ versus \z\ of pulsars, which have either measured parallaxes or are associated with globular 
clusters at known distances. DM ~ n e l is the dispersion measure of a pulsar at a distance I. 

The NE2001 model w as widely used, although shortcomings wer e noted so on such as the d e crease 
of z with pulsar distance dKramer et all 120031: lLorimer et al.ll2006l) . Recently, iGaensler et"aD (l2008h 
convincingly demonstrated that many DMs of low-latitude pulsars are influenced by local structures 
such as H II regions along the line-of-sight and derived a scale-height of about 1.8 kpc by inclu ding 
only pulsars at \b\ > 40°. The mid-plane density derived was 0.014 cm ~ 3 . IGaensler et al.l (l2008h also 
found that the filling factor of the thermal gas, which was not addressed in the NE2001 model, grows 
exponentially in the rang e of < \z\ < 1.4 kpc w ith a scale-height of 0.7 kpc. Th is differs from the 
filling factors obtained bv lBerkhuijsen et al l (l2006l) and lBerkhuij sen & Mullerl(l2008l) . 
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Measured distance (kpc) Measured distance (kpc) 

Fig. 1 Measured distances versus distance estimates from the NE2001 model (left panel) 
and the modified NE2001 model (right). The measurements from parallaxes are indicated by 
open circles and those from the associated globular clusters by filled circles. The lines mark 
the case that the measured and estimated distances are equal. 

In the following we modify the NE2001 model by re placing the sca l e-heig ht and the mid-plane 
density for the thick disk component by those provided by IGaensler et al.l d2008l) . We investigated the 
modified NE2001 model by comparing the distances of pulsars estimated from the model by their quoted 
DMs with that from independent measurements such as parallax a nd the associated glo bular clusters. 
The up-to-date parallax measurements for pulsars were collected bv lVerbiestet"aHd2010l) . We obtained 
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DMs of these pulsars from the ATNF pulsar database^. Data for pulsars in globular clusters were taken 
from the web-pag^l maintained by Paulo Freire, where the distances of the clusters and the DMs of 
pulsars are listed. The results from the original as well as from the modified NE2001 model are shown 
in Fig.Q] Although there is a considerable scatter, the distances of pulsars estimated from the modified 
NE2001 model agree better with the observations than the original NE2001 model, especially for pulsars 
with distances larger than about 6 kpc. 

We do not claim that our modified NE2001 model is a complete substitute for the NE2001 model 
in describing the diffuse ionized gas in the Galaxy. Other NE2001 components such as the thin disk 
and the arm parameters need to be revised correspondingly to interpret properly interstellar scattering 
and scintillation. A new model, NE2008, accounting for all new relevant observation is currently been 
developed (Jim Cordes, private communication). 
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Fig. 2 Longitude and latitude profiles from the simulated and the template free-free emission 
maps at 22.8 GHz. The thin solid lines are from the WMAP five-years ME M template, the 
thick s olid lines from the modified NE2001 model with the filling factors bv lGaensler et al.l 
(120081). and the dotted li nes from the original NE2001 model with filling factors derived by 
iBerkhurjsen et al.l d2006l) . 



W ith the modified NE2001 thermal electron density model and the filling factors bv lGaensler et al.l 
d2008l) as i nput w e simulated the all-sky free-free emission at 22.8 GHz. This i s compared to the result b y 
I Sun et al" I (120081) based on the original NE2001 model and t he filling factors b y lBerkhuiisen et al I (120061) . 
The WMAP five-years MEM thermal emission template dGold et all 120091) is shown for comparison. 
Slices extracted from the above mentioned maps are shown in Fig. [2] Note that the lower envelope of the 
observed emission is to be compared with the simulated profiles, which do not include all individual H II 
regions nor other local featu res. Along the Ga lactic plane thermal emission from the modified NE2001 
model is lower than that bv lSun et al.l d2008l) . but still agrees with the observations. The reduction of 
thermal emission results from the smaller mid-plane electron density of 0.014 cm~ 3 in the modified 
model and 0.034 cm" 3 in the NE2001 model. 



2 http://www.atnf.csiro.au/research/pulsar/psrcat/ 

3 http://www.naic.edu/~pfreire/GCpsr.html 
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The revised scale height by IGaensler et al.1 ( 2008 ) was not accepted everywhere. I Savage & Wakkerl 



(120091) used the same data as IGaensler et al.l d2008l) . but invoked a different fitting scheme, where 
an adjustable pa tchiness error was includ ed to assure that xl ~ 1- The fit by IGaensler etail (l2008h 
yie lded xl ~ 5.ISavage & Wakkerl d2009l) obtained a scale-height of about 1.4 kpc instead of 1.8 kpc 
by IGaensler et al.l ( 2008 ). Their corresponding mid-plane density is then about 0.016 cm~ 3 . We 
run simulations by mod ifying the thick disk components of the NE2001 model with the values by 
ISavage & Wakkerl (120091) . and found the results do not vary significantly from those with a scale-height 
of 1.8 kpc, which we presented above. New pulsar data became recently available being almost con- 
sistent with the 1.8 kpc scale-height also when applying the fitting method used by ISavage & Wakkerl 
(120091) (Bryan Gaensler, private communication). 



3.2 Regular magnetic field properties 



It is a customary to split the regular Galactic magnetic field into a disk and a halo component. The 
configuration of the disk fields were usually classified into three types, (1) axi-symmetric spiral (ASS); 
(2) bi-symmetric spiral (BSS); (3) following the spiral arms. The commonly used halo field patterns are 
toroidal, poloidal and a combination of both. The disk and halo fields and in particular their direction 
were mainly constrained by observed RMs. 



3.2.1 The disk field 



Early models for the disk magnetic field frequently used to study the propag ation of u l trahig h energy 
cosmic- rays do not agree wit h the RMs observ ed from extragalactic sources dSun et all 120081) or from 
pulsars dNoutsos et al.Ll2008l) . ISun et al.1 d2008l) have proposed three new models for the disk magnetic 
field configuration: ASS+RING with reversals in the Galactocentric rings, BSS, and ASS+ARM with 
reversals in arms. All models predict RMs consistent with the available data. Reversals of the large-scale 
magne tic field have been subject of a controversial debate for a long time. Recently, iNota & Katgertl 
(120101) concluded that the magnetic field in the fourth Galactic qua drant exhibits rev ersals at the arm- 
interarm interfaces, which confirmed the earlier results obtained bv lHan et al.l d2006l) . However, given 
the uncertain distances of pulsars and the difficulties to properly assess the influence of foreground 
structures such as H II regions, supernova remnants, Faraday Screens or the giant local loops, this topic 
clearly needs further investigations. 

Here we kept the disk field of the SRWE08 models unchanged and simulated the RM map with 
the ASS+RING model as an example, which was recently further supported by new RM measurements 
(Ivan Eck & Browr3 , l2009l) . The RM profile along the Galactic plane taken from the maps is compared to 
the binned CGPS and SGPS RM data (Fig. [3] upper-left panel). The RM profile based on the modified 
NE2001 model and the ASS+RING disk field basically reproduce all the observ ed features . In gen eral, 
the RM magnitudes from the present model are smaller than those modelled by ISun et aL (2008). For 
example, the RM gradient in the range of 80° < I < 150° is now more shallow than before. To reach 
an agreement between the updated and the SRWE08 models we need to increase either the disk mag- 
netic field strength or the thermal elect ron density . The l o cal field strength o f about 2 fj,G is determined 
by the RM/DM ratio of pulsars (e.g. lHan et al 1 120061) . iBeck et al.l d2003l) pointed out that the esti- 
mate of magnetic field strength by the RM/DM ratio is biased if t here exists a corr elation between the 
thermal electron density and the magnetic field strength. However. IWu et all (120091) performed simula- 
tions which showed that the RM/DM ratio correctly represents the mean strength of the magnetic field. 
Although there remain uncertainties on the turbulent properties of the interstellar medium, we consider 
the 2 /iG local field strength as very robust. By increasing the mid-plane thermal electron density by 
about a factor of two it is about the same as in the NE2001 model and the RM profile from the present 
model is consistent with that from the SRWE08 model. However, the DM sin \ b\ versus z could not be 
fitted well with the larger mid-plane electron density and a scale-height of 1.8 kpc. Possibly the ex- 
ponential description of the thermal electron density is oversimplified or an additional electron density 
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Fig. 3 RM profiles along Galactic longitudes and latitudes. The solid lines are from the 
present ASS+RING model, and the dashed lines from the correspondi ng SRWE08 model . 
The R Ms in the plane (upper-left panel) are from the CGPS and SGPS dBrown et al.ll2003l 
120071) . RMs in the area 100° < I < 120° and -5° < b < 20° (upper-right panel) are from 
the CGPS extension (Jo- Anne Brown, in prep.), and in the halo regions from NVSS sources 
dTavloret all 120091) . 



component exists, which contributes just near to the plane. We expect this point will be solved by the 
expected NE2008 model. 

3.2.2 The halo field 

The asymmetric distribution of RMs in longitude and latitude relative to the Galactic plane and the 
Galactic centre indicates that the Galactic halo field has opposite signs below and above the plane. 
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ISun et al.l (120081) have described the toroidal halo field B<p(R,z) following lProuza & Smfdal (120031) by, 

1 R f R — Rq 

\z\ — z \ K Q \ R 



B<p(R,z) = sign(z) J B — — — cxp( - ), (1) 



1 



Zl 



where sign(z) takes the signs of z. 

To account for the revised thermal electron density model, we obtained new parameters for the 
halo field as: z = 1.5 kpc, Z\ = 0.2 kpc for \z\ < z and otherwise Z\ = 4 kpc, B = 2 fiG, 
and Rq = 4 kpc. From new all-sky simulations with added halo and disk magnetic fields, we obtained 
the RM latitude profile for the longitude range of 100° < I < 120° (Fig. [3] the upper-right panel) to 
be compared with RM data from the CGPS high-latitude extension (Brown et al., in prep.). Near the 
Galactic plane the simulated RM profile does not agree with the observations, which again is caused by 
the too small mid-plane thermal electron density in our modified NE2001 model. The deviations vanish 
by increasing the thermal electron density by a factor 2 as mentioned in Sect. l3.2.Tl to reconcile the RM 
difference in the plane. 

We also obtained average RM longitude profiles and show them for 25° < b < 35° and 40° < 
b < 50°, as well as their southern counterparts i n Fig. [3] (middle an d lower panels). These profiles were 
compared with the RMs from NVSS sources dTavlor et all l2009h binned in 10° longitude intervals. 
RM profiles from new simulations generally reproduce the observed RM anti-symmetry with smaller 
amplitudes compared to the SRWE08 models. 

There is evidence based on polarization data from the Global Magneto-Ionic Medium Surve y 
(GMIMSfl currently running at the DRAO 26-m telescope at L-band dWolleben et all 120091 1201 Obi) , 
that magnetized filaments are associated with a very local and therefore very extended H I shell located 
above the Galactic plane in the direction of the Galactic centre w ith large positive and negative RMs in 
the first and forth quadrant, respectively dWolleben et all 1201 Oal) . Very likely this bubble explains the 
RM deviations between the current halo model and the RM measurements seen in Fig. [3] for the area 
30° < I < 60° and 25° < b < 35°. This finding does not question the present halo model, as its 
asymmetry above the plane of 10-20 rad m -2 is smaller t han that attributed to the magnetized filaments. 
In the region 45° < I < 60° and -60° < b < -20° iTavlor et alj d2009l) derived positive RMs for 
distinct small areas (see their Fig. 4), which partly cancel the in general negative RMs in this area and 
cause deviations between simulations and observations. These small areas with positive RMs may be 
associated to southern Loop I filaments. 

It is important to keep in mind that we observe the halo magnetic field in superposition with the 
disk field. Because the toroidal magnetic field above the plane is directed opposite to the disk field, 
the intrinsic RM asymmetry from the halo model is reduced. Vice versa, below the plane disk and 
halo field add and the halo asymmetry becomes more evident (Fig. [3]). It is therefore very important 
to obtain a large set of RM data for the southern sky, especially for the region of 240° < I < 360°, 
to prove the predicted halo RM asymmetry. We note that the northern sky area with large negative 
RMs within 60° < I < 140° below the Galactic plane towards b ~ —40°, named region A by 
ISimard-Normandin & Kronberd dl980l) . might be influenced by Loop II, which is generally believed 
to b e a local shell struc t ure. 

iNagar & Matulichl d2010h suggested that a vertical field component needs to be added to the 
SRWE08 models to better describe deflections of ultrahigh energy cosmic-rays. This vertical field 
component has a strength of about 0.2 fiG at the solar position consistent with that derived by 
lHan & Oiaol d 19941). We p erformed simulations including a vertical field of that strength as proposed 



by Prouza & Smfdal d2003l) and found the results just marginally deviate from the simulations without 



this com ponent. Thi s mean s that such a weak vertical field cannot be constrained by our simulations. 
Recently, iMao et "aH d2010h obtained RMs of numerous extragalactic sources towards the northern and 
southern Galactic pole and convincingly demonstrated that there does not exist a coherent vertical field 
at the Sun's position. Therefore we do not include any vertical magnetic field component in our models. 



https://www.as trosci.caAisers/drao/gmims 
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3.3 Random magnetic fields and the cosmic-ray electron density 

3.3.1 The random magnetic field component 

I Sun et all (120081) and ISun & Reichl (120091) have invoked isotropic random magnetic fields with a 
Gaussian distribution to simulate all-sky maps. Sim ulations of high angular resolution patches of the 
sky follow a power-law spectrum. Ijaffe et all (1201 d) have included an additional magnetic field com- 
ponent, called "ordered" component, which is a regular field with numerous small-scale field reversals. 
Unlike isotropic random fields, this component contributes to total and polarized emission, which in 
turn reduces the cosmic -ray electron density to match observations. The "ordered" magnetic field does 
not increase RM, but increases its scatter. Although this "ordered" component is clearly helpful to reach 
consistency between simulations and observations, the spatial scales of its reversals are entirely unclear 
and there is no proof so far on its existence in our Galaxy at all. We therefore do not include an "ordered" 
component in our simulations. 




b (=) b (=) 



Fig. 4 Total intensity profiles at 408 MHz. The thin solid lines s how profiles obtained from 
the 408 MHz all-sky total intensity survey (lHaslam et all Il982l) . The results from the new 
model and the SRWE08 model are displayed as thick solid and dashed lines, respectively. 
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3.3.2 The cosmic-ray electron density distribution 

The SRWE08 models truncated the distribution of the cosmic -ray electron density at \z\ = 1 kpc to 
avoid excessive synchrotron emission at high latitudes caused by the strong halo magnetic field. This 
truncation is quite certainly physically unrealistic. Since the halo magnetic field in the revised models 
is much smaller thi s cutoff in z is not needed any more. A scale-height of 0.8 kpc instead of 1 kpc used 
bv lSun et ai] (l2008[) was found to adapt better to the data. 

With the revised models, we simulated an all-sky total intensity map at 408 MHz and a polarized 
intensity map at 22.8 GHz at an angular resolution of 15'. The 4 08 MHz map was the n smoothed to 51', 
which is the angular resolution of the 408 MHz a ll-sky survey by Haslam et al.l d 19821) . The simulated as 

20091) " were smoothed to 2° to increase 



well as the WMAP 22.8 GHz polarization maps dHinshaw et al 



the signal-to-noise ratio. Slices extracted from the simulated maps are shown in Fig. [4] for 408 MHz 
total intensities and in Fig.[5]for 22.8 GHz polarized intensities. The new simulations qualitatively agree 
with the observations at a level comparable to the SRWE08 models. 
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3.3.3 The local enhancement 

The SRWE08 models include a local emission excess, realized by enhanced cosmic -ray electrons within 
1 kpc, to account for the increase of the synchrotron emissivity towards the solar system based on low- 
frequency absorption data from optically thick HII regions. Clearly the available data are quite limited 
and not suited to map the local synchrotron excess in 3D. The spherical approach used in the SRWE08 
models is therefore almost arbitrary. It will be the task of LOFAR or other new low-frequency telescopes 
under construction to observe many more faint H II regions across the sky to model the local synchrotron 
excess in 3D. This task is closely related to obtain a realistic halo model as both components are present 
at high latitudes. We note from Fig. g] that for the region of 180° < I < 270° the observed 408 MHz 
total intensity at high latitudes is below that expected from the simulations. This indicates a possible 
offset of the local excess centre from the solar position. We shifted the centre of the local enhancement 
towards I = 45°, b = 0° with a distance of about 560 pc from the Sun and repeated our simulations. The 
results are shown in Fig. [6] Quite obviously a much better agreement between the simulation and the 
observations is achieved shifting the centre of the local enhancement. However, there are not yet clues 
from observations on the centre location and the 3D geometry of the local enhancements, so that we are 
looking forward to numerous H II region absorption results obtained by LOFAR. 




-90 -60 -30 30 60 90 -90 -60 -30 30 60 90 



b (°) b (°) 



Fig. 6 Total intensity latitude profiles. The thin lines are from the 408 MHz all-sky survey. 
The solid lines are from the simulations with the local enhancement shifting about 560 pc 
from the Sun, and the dashed lines are from the SRWE08 model. 



4 SUMMARY 



In this paper we updated the parameters of the 3D-emissivity models used bvlSun et alj d2008l) fo r all- 



sky simulations based on a revised thermal electron density scale-height bv lGaensler et al.l (|2008i). The 



scale-height and mid-plane electron density of the NE2001 thermal electron density model of 1 kpc and 
0.034 cm -3 were replaced by 1.8 kpc and 0.014 cm -3 . In consequence, the maximum halo field strength 
of 2 /iG, and a scale-height of the cosmic -ray electron density of 0.8 kpc is obtained, which we consider 
to be much more physically relevant than the 10 ^G field strength and the truncation with z at 1 kpc 
obtained by lSun et al.l d2008l) . With these modifications, we reproduce the high latitude RM distribution 
of extragalactic sou rces, total intensity and polarized intensity all-sky survey maps in the same way 
as ISun et all d2008l) . We note that the mid-plane thermal electron density needs to be increased by a 
factor of about two to properly represent the RMs in the Galactic plane. We also find that a shift of the 
centre of the local synchrotron enhancement from the Sun position by 560 pc towards I = 45°, 6 = 0° 
significantly improves the models to adapt to total intensity observations towards the outer Galaxy. 
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